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FREE JET IMPINGEMENT NORMAL TO A CURVED SURFACE IN A VACUUM* 

By Leonard V. Clark 
Langley Research Center 

SUMMARY 

An experimental investigation of highly underexpanded exhaust jets impinging nor- 
mal to a curved surface was conducted in the Langley 55-foot vacuum cylinder. The 
investigation included measurements of impingement pressures and of various bow- 
shock characteristics. The curved surfaces were semicircular cylinders and impinge- 
ment pressures were measured both along and around the cylinders which were located 
at various distances from the exhaust jets. Three nozzles were used with unheated 
nitrogen and helium; thus, six different sonic and hypersonic jets could be studied. The 
tests were conducted at pressure altitudes at or above 200,000 feet (61 km). The exper- 
imental data are compared with available theory to determine the applicability of analyt- 
ical methods for predicting exhaust -jet impingement pressures for a range of conditions. 

The results of this investigation verify previous experience and analysis and indi- 
cate that Roberts' simplified analysis is as good as the more cumbersome method of 
characteristics plus Newtonian analysis for determining impingement pressures resulting 
from underexpanded exhaust jets. 


INTRODUCTION 

Exhaust jets are finding many useful applications in space at the present time such 
as attitude control of spacecraft and spacecraft propulsion. The expansion of these jets 
because of the low ambient pressure in space may complicate such operations as staging, 
rendezvous, and docking, and even a lunar landing. During orbital rendezvous of vehicles 
with space stations, for example, the vehicle will necessarily maneuver in many direc- 
tions. The possibility of the exhaust from thrust or control systems being directed 
toward the station is not unlikely. If the exhaust impinges upon the station, it may pro- 
vide a sufficient force to change the stations' orbital characteristics. During a lunar 

*The information presented herein is based upon a thesis entitled "An Investigation 
of Highly Underexpanded Exhaust Jets Impinging Normal to a Curved Surface" submitted 
by Leonard V. Clark in partial fulfillment of the requirements for the degree of Master of 
Science in Engineering Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, 
March 1966. 
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landing, a vehicle like the Apollo lunar module (LM) must brake its approach velocity by 
using retrograde rockets. One of the problems to be considered in this instance is the 
exhaust plume interference with the landing radar system as the LM maneuvers in closing 
with the lunar surface. These and other problems are summarized in reference 1. 

The characteristics of exhaust jets depend to some extent upon the ratio of jet-exit 
pressure to ambient pressure. When the ambient pressure is less than the jet-exit pres- 
sure, some expansion takes place outside the jet nozzle. A jet operating under this con- 
dition is generally referred to as an underexpanded jet. The character of underexpanded 
free jets is well known as evidenced by references 2 and 3, two of many published works. 
However, a knowledge of the character of impinging jets has only recently become impor- 
tant and is less well known. There are few published experimental papers and the ade- 
quacy of available theory for predicting impingement pressures has not been widely 
demonstrated. The studies of references 4 to 6 were concerned with experimentally 
determining impingement pressures on flat plates oriented parallel to or at slight angles 
to the jet exhaust axis. The applicability of an analysis using the method of characteris- 
tics and Newtonian flow theory (hereinafter called characteristics plus Newtonian analy- 
sis) was demonstrated in reference 5 for jets at one location parallel to a plate and using 
nitrogen and helium. The studies of references 7 and 8 were concerned with experi- 
mentally determining impingement pressures on flat plates oriented normal to the jet- 
exhaust axis. Variables, such as distance of the jet away from the plate, jet-exit Mach 
number, nozzle contour, and the ratio of jet-exit pressure to ambient pressure, were 
studied. The applicability of a characteristics plus Newtonian analysis and an approx- 
imate analysis developed by Roberts in references 9 and 10 was demonstrated in 
reference 8. 

An examination of available literature revealed little information regarding the 
influence of either scale effect of jet size relative to the surface or surface shape on the 
distribution of impingement pressures. The present study was undertaken to determine 
the applicability of available theory for predicting impingement pressures for a range of 
these conditions. Impingement pressures were measured on several semicircular cyl- 
inders having different curvature for several under expanded sonic and hypersonic exhaust 
jets. The jet axis was maintained perpendicular to the surface longitudinal axis of sym- 
metry to provide data directly comparable with theory. The experiments were conducted 
in the Langley 55-foot vacuum cylinder and included a shadowgraph study of the exhaust- 
impingement region on one of the surfaces for the different jets from which certain bow- 
shock characteristics were measured. 
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SYMBOLS 


The units used for the physical quantities defined in this paper are given both in the 
U.S. Customary Units and in the International System of Units (SI). Factors relating the 
two systems are given in reference 11. 

A area 

d diameter 

dA differential area 

F thrust 

h distance measured along jet axis from jet exit to surface 

K constant given by equation (11) 

l distance measured along model surface in axial direction from jet-axis 

impingement point 

M Mach number 

p static pressure 

q dynamic pre ssure 

r radial distance measured from jet axis 

R gas constant 

s distance measured along model surface in circumferential direction from jet- 

axis impingement point 

T temperature 

V velocity 

x distance measured along jet axis from jet exit 

a streamline inclination angle measured from jet axis 
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y ratio of specific heats 

6 initial expansion of exhaust at nozzle exit measured between jet axis and 

tangent to jet boundary 

A bow -shock standoff distance along jet axis 

6 azimuthal angle measured from jet center line at jet exit 

d Q nozzle -exit half -angle 

v Prandtl -Meyer expansion angle 

p density of exhaust gas 

a> angle measured from jet axis at model center of curvature 

Subscripts: 
a jet axis 

B bow -shock bowl 

c nozzle settling chamber 

e nozzle exit 

m model 

max maximum 

N normal 

r shock recovery conditions 

s surface 

stag jet-axis impingement point (stagnation) 

T tangential 
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1 denotes a particular point in the exhaust flow field 

00 vacuum chamber conditions (ambient) 

APPARATUS AND PROCEDURE 
Test Setup and Procedure 

The experimental investigation was conducted in the Langley 55-foot (16.7 m) 
vacuum cylinder with the test setup shown in figure 1. High-pressure unheated gas was 
supplied to the nozzle from a bottle farm located outside the cylinder with the nozzle - 
chamber pressure not regulated. Nozzle -chamber pressures measured during the 
experiments ranged from 340 to 1630 psia (2.34 x 10® to 11.2 x 10® N/m 2 ) and the 
vacuum -chamber pressure prior to initiation of jet-exhaust flow was between 0.050 to 
0.155 torr (6.67 to 20.6 N/m 2 ). These nozzle -chamber pressures are realistic of those 
for present launch -vehicle propulsion systems while nozzle -chamber pressures of 
approximately 100 psia (0.689 x 10® N/m 2 ) are realistic for operation of exhaust jets 
for spacecraft control and propulsion. The vacuum -chamber pressure would correspond 
to an altitude above the earth of approximately 200,000 feet (61 km). Figure 1 shows the 
general arrangement of the nozzle, quick-acting valve, impingement surface, shadow- 
graph light source, mirrors, and 70-mm sequence camera. 

The impingement surfaces, semicircular cylinders, were mounted perpendicular 
to the nozzle axial center line. Surface pressure distributions were measured both 
along and around the cylinders which were located at various distances from the six dif- 
ferent exhaust jets used in this investigation. A solenoid valve located in the gas supply 
line a short distance from the nozzle permitted the nozzle to be operated remotely. In 
addition, the test nozzle was mounted on a movable carriage; this mounting permitted the 
distance between the nozzle and the surface to be varied remotely. These features made 
vacuum -chamber repressurization necessary only for model or nozzle changes. The 
duration of a typical test at a particular jet height was about 3 seconds during which time 
the vacuum chamber pressure increased by about 0.150 torr (20.0 N/m 2 ). For a sub- 
sequent test, the nozzle height was changed and the volume of gas exhausted into the 
vacuum chamber from the previous test was evacuated by pumps which were kept in con- 
tinuous operation. 


Impingement Surfaces 

The cylindrical impingement surfaces are shown dimensionally in figure 2. The 
locations of the static -pressure orifices along and around the surface of the cylinders 
are indicated in figure 2. Also shown is the coordinate system l , go at the center of 
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which the jet-exhaust axis was alined. The static -pres sure orifices were 0.038 inch 
(0.0965 cm) in diameter. 


Test Nozzles 

The investigation was conducted with three nozzles which are shown dimensionally 
in figure 3. Using two gases with each of the nozzles actually permitted six different 
exhaust jets to be studied. One of the nozzles was a converging (M e = 1.0) nozzle 
whereas the other two nozzles were converging-diverging conical nozzles with an area 
ratio of 25 with an exit half-angle of 15°. The two converging-diverging nozzles differed 
only in size, one being twice as large as the other. The nominal design exit Mach num- 
ber of the converging-diverging nozzles was 5.0 with nitrogen and 7.09 with helium. The 
investigation reported in reference 6 used a similar converging-diverging nozzle with air 
and showed by a static -pressure measurement in the expansion wall near the nozzle exit 
that the jet actually had an indicated Mach number of 4.79 instead of the inviscid design 
value of 5. In addition, measurements of the initial turning angle of the exhaust flow at 
the nozzle exit 5 obtained from schlieren photographs at known values of the ratio of 
nozzle chamber pressure to ambient pressure indicated an effective nozzle-exit half- 
angle 5 of 26.5° instead of the inviscid design value of 15°. 

The calculated exhaust initial turning angle for the jets used in this study for var- 
ious values of the ratio of jet -exit pressure to ambient pressure is shown in figure 4. 

This angle was computed by using the approach discussed in reference 1 for the prop- 
erties of nitrogen and helium ascertained from the compressible flow tables of refer- 
ences 12 and 13. Also shown in the figure is the maximum possible turning angle for 
each of the nozzles which occurs at an infinite pressure ratio. The figure indicates that 
even at the finite operating pressure ratios of the present tests, the jets were not too far 
from being fully expanded. 


Instrumentation 

Nozzle -chamber pressure was measured by a 2000-psia (13.8 x 10® N/m^) capacity 
pressure transducer located in the nozzle -settling chamber as shown in figure 3. The 
vacuum -chamber ambient pressure was measured before and after the short test period 
by a thermopile -type pressure gage that had been calibrated with a McLeod gage. The 
impingement surface was equipped with static -pressure orifices which were connected to 
strain-gage diaphragm -type absolute pressure gages by short lengths of plastic tubing. 

The orifice on the jet axis was connected to a 10-psia (68.9 x 10® N/m®) capacity gage. 
The surrounding four orifices were connected to 5-psia (34.5 x 10® N/m^) capacity gages. 
The next orifice in each row was connected to a 3-psia (20.7 x 10® N/m^) capacity gage 
and the remaining orifices were connected to 1-psia (6.89 x 10® N/m2) capacity gages. 
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All pressures were recorded during the tests on an oscillograph. Nozzle position was 
indicated on the same oscillograph by the output of a potentiometer. During the tests the 
pressure gages were operating anywhere from 10 to 100 percent of the gage rating and 
it is estimated that the overall accuracy of the data is ±2 percent. 

The response time of the various pressure gages is of importance because of the 
short test time. In figure 5 are presented the nozzle -chamber and surface stagnation 
pressure history for a portion of a typical test. For the present tests the impingement 
pressure data were selected at a time when a constant ratio of surface to nozzle -chamber 
pressure existed. This time ranged from 0.3 second for the gage measuring stagnation 
pressure (fig. 5) to 0.7 second for the outermost gages. 

During some of the tests, shadowgraphs of the jet exhaust in the vicinity of the 
impingement surface were obtained by using the system depicted in figure 1. For this 
system a beam of light produced by a commercial, 25-watt zirconium concentrated-arc, 
point source located at the focus of an 18 -inch -diameter (0.457 m) parabolic mirror was 
adjusted to be parallel to the impingement surface. This beam was intercepted by 
another 18 -inch -diameter (0.457 m) parabolic mirror located on the other side of the 
vacuum chamber and was directed into a lensless 70-mm sequence camera. Both the 
light source and camera were contained in pressure boxes which were vented to atmos- 
pheric pressure and equipped with glass viewing ports. 

ANALYSIS 

The experimental data obtained during the present tests were compared with the 
results obtained from two methods of approximate analysis. Both methods use Newtonian 
flow theory to get the impingement pressure; however, the properties of the jet exhaust 
are obtained differently. One uses the method of characteristics whereas the other uses 
a simplified integral method. A good discussion of Newtonian flow theory is given in 
reference 14. 


Characteristics Plus Newtonian Analysis 

In this method the flow properties of the exhaust are calculated by the method of 
characteristics, and then the surface pressures are estimated by Newtonian theory. A 
discussion of the method of characteristics may be found in reference 15 which also 
provides a computer program for calculating the characteristic network for any jet 
operating in an underexpanded condition. This particular program yields jet-exhaust 
boundary coordinates and several additional internal parameters in the flow field, one of 
which is the local Mach number. To determine impingement pressures, plots of constant 
Mach number contours and constant-flow-inclination contours within the exhaust plume 
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were required for use with Newtonian theory. A computer program to interpolate 
between consecutive points in the characteristic net to obtain these contours was avail- 
able only for Mach number contours (ref. 16). Available calculated contours were 
utilized, where possible, for expediency. For example, reference 8 contains the Mach 
number contours for a sonic jet and an M g = 4.79, = 26.5° jet for air (y = 1.4). 

These contours, shown in figure 6, were considered to be applicable for the nitrogen 
jets of the present study as the properties of nitrogen and air (79 -percent nitrogen) 
are very similar. They were computed for a specific ratio of nozzle -chamber pressure 
to ambient pressure; however, as noted in reference 8 the contours are unaffected by 
additional increases in pressure ratio. In figure 6(b) can be seen the leading character- 
istic line as determined from the nozzle -exit Mach number by the method detailed in 
reference 2. The leading characteristic line serves as a dividing line between the 
internal flow region where the flow characteristics are determined by the nozzle contour 
and the external flow region where the flow characteristics are determined by the ratio 
of jet -exit pressure to ambient pressure. For conical nozzles the leading characteristic 
line is curved and the flow in the internal region is radial and axisymmetric. The con- 
stant Mach number contours for an M g = 7.09, = 15° helium jet were obtained by 

using the computer programs in references 15 and 16 and are given in figure 7. The 
flow-inclination contours for the nitrogen jets shown in figure 8 were calculated for but 
are not presented in reference 8. These contours are seen to be approximately equal to 
azimuthal angles measured from the jet axis at the jet exit. 

Values of p, p r , and q at any point in the jet -exhaust flow field may now be 
found from the isentropic relations tabled, as a function of Mach number, in references 12 
and 13 for air and helium after which surface pressures are obtained by using Newtonian 
flow theory, a shock layer coincident with the surface being assumed. By utilizing the 
nomenclature shown in sketch (1) which depicts a typical jet exhaust streamline Vj 
impinging upon a body with a circular cross section, 



Sketch (1) 
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the impingement pressure is obtained from Newtonian theory as 


p s 



+ u> 



( 1 ) 


where 


ii = | p i v i 2 < 2 > 

Along the model longitudinal axis, the geometry simplifies the impingement pressure 
equation to 


P s = 2q x (cos a^ 2 + p x 

At the jet axis the impingement pressure equation would further simplify to 


(3) 


p stag = 2q a + p a (4 ' 

The static pressure p^ may generally be neglected when the Mach number of the inci- 
dent flow is very great as this term is much smaller than q^ for large M^. This 
statement means that the actions exerted by the fluid particles upon each other in the 
impinging flow are small in comparison with the actions exerted by the same particles 
on the body surface when the flow is highly supersonic. Even the flow from a sonic jet 
becomes highly supersonic downstream of the nozzle exit as it expands to the vacuum; 
thus the surface pressure at the jet axis may be written 


p stag ^a 


(“» » l) (5) 


Another approach based upon Newtonian flow theory is given in reference 17 where 
it is assumed that the surface pressure at a point is equal to the bow-shock recovery 
pressure for the local flow Mach number at that point multiplied by the cosine squared of 
the streamline impingement angle. Thus, for the curved surfaces, the impingement 
pressure may be written as 


p s = p r,l[ cos (“l + w l)] 

The surface pressure along the model axis then becomes 

p s = p r,l( COSO! l) 2 


( 6 ) 

(7) 
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and the surface stagnation pressure is given by 


p stag p r,a ^ 

To apply Newtonian theory along the jet axis requires only a plot of the axial vari- 
ation of the flow-field Mach number for the particular jet. Figure 9 shows the calculated 
axial variation of Mach number for some of the present jets. Included is the variation 
for an M g = 5, 6 & = 15° air jet which was obtained from reference 8. A comparison 
of jet-axis impingement pressures obtained from the two approaches using a character- 
istics plus Newtonian flow analysis (eq. (5) and eq. (8)) is afforded by figure 10. In this 
figure calculated values are plotted for various values of Mach number for air and 
helium. The figure shows that surface pressure is expected to decrease with increasing 
flow Mach number and that both approaches give about the same results. 

Analysis of Roberts 

In this method the flow properties of the exhaust are calculated by an approximate 
integral method, and then the surface pressures are estimated by Newtonian theory. This 
method is developed in references 9 and 10 and considers two distinct cases of flow: 

(1) flow where the leading characteristic intersects the jet axis somewhere downstream 
(called typical exhaust flow) and (2) flow where the leading characteristic does not inter- 
sect the jet axis (called spherical source flow). The two flow cases are distinguished 
according to reference 9 by whether or not 


v + 29 2 v =90 
e' e max 



( 9 ) 


This equation is graphically represented in figure 11 and delineates the maximum jet- 
exit Mach number for which the leading characteristic of a given conical nozzle can inter- 
sect the jet axis. 


Typical exhaust flow .- For this flow case the governing impingement pressure 
equation when the jet is far from the curved surface is given by the method of refer- 
ence 10 as 


where 



K+ 2 




( cos 0) K+2 [cos(0 + «)] 
p c A e ~ ~ 


K = y(y - l)M e 2 


( 10 ) 


(ID 
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F is the jet thrust, and 6 is the azimuthal angle measured from the jet axis which 

locates a particular point in the flow field. This equation has the restriction that 
2 

M g » 1; and, in addition, the assumption was made that the exhaust gas velocity is con- 
stant and equal to V max where 


V 


max 



( 12 ) 


which is not exactly true near the jet exit. Along the model longitudinal axis of symme- 
try, the impingement pressure equation would simplify to 


K + 2 / h 
2 r. 


F <cos 0) K+4 


PrA 


(13) 


"c " V e / ~e 
At the jet axis, the impingement pressure equation would further simplify to 


p stag = K + 2 / h \~ 2 F 

p c 2 \ r e J P c A e 


(14) 


When the jet is close to the impingement surface, reference 9 gives the surface 
pressure as 


p stag _ F _ p e^e + p e 
p c P C A e Pc 


(15) 


The equations for impingement pressures for the different jet positions relative to the 
surface (eqs. (14) and (15)) are seen to coincide at 


h 


r 


e 


K + 2 


(16) 


which is taken as the limit of applicability of equation (14). 

Spherical flow .- For this flow case reference 9 gives an asymptotic equation for 
the density distribution of the exhaust gas along the jet axis but does not present any 
equations for determining impingement pressures. By removing the restriction on the 

jet height ratio that — » 1 and using Newtonian theory with this density distribution, 

r e 

the equation which gives surface stagnation pressure for spherical flow jets becomes 
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The other assumption made during the derivation of the equation of axial density distri- 
bution for these jets was that the exhaust gas velocity is constant and equal to V mQv . 

RESULTS AND DISCUSSION 

Impingement Pressure Data 

Experimentally determined pressure distributions both around and along the various 
curved models resulting from the different impinging jets of the present study are given 
in figures 3.2 to 19 for several representative jet heights. Figures 12 to 15 show pres- 
sure distributions around the cylinders and figures 16 to 19 show the associated pressure 
distributions along the cylinders. 

Effect of Jet Distance From Surface 

The data of figures 12 to 15 show that the distribution of impingement pressures is 
generally symmetric about the jet-axis impingement point and that surface pressures 
increase as the jet approaches the surface. There was some nozzle misalinement for 
one configuration (data of fig. 12(a)) which was corrected for subsequent tests. The var- 
iation of the axial impingement pressure with jet distance from the surface, which is 
generally the maximum surface pressure, as obtained from figures 12 to 15 is summa- 
rized in figures 20 to 24. 

M = 5 nitrogen jet.- Figure 20 is a summary of the experimental data for both 

V 

sizes of the M 0 = 5 nitrogen jet and includes analytical estimates as well as experi- 
mental data from references 7 and 8 for the same Mach number jet by using air impinging 
a flat plate. Axial impingement pressure depends on the ratio of the distance of the jet 
from the surface to the jet-exit size (jet height ratio) since the data from three different 
sized jets fall along a common path. Axial impingement pressure does not depend upon 
the shape of the surface since the data from four different curved surfaces and aflat 
plate fall along the same common path. Analytical estimates from characteristics plus 
Newtonian analysis (eqs. (5) and (8)) utilizing the appropriate axial Mach number distri- 
bution from figure 9 agree very well with the experimental data especially when the jets 
are close to the surface. This analysis slightly underestimates the decrease of surface 
stagnation pressure for jets at large distances from the surface. For the analysis based 
on Roberts' approach, the equations for a typical exhaust flow case are applicable. (See 
fig. 11.) When the jet is within a couple of jet-exit diameters from the surface, the axial 
impingement pressure is given by equation (15) and at further distances from the surface, 
the axial impingement pressure is given by equation (14). This analysis agrees very 
well with the experimental data when the jet is 4 to 10 exit diameters from the surface. 
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Noted on this and subsequent similar figures is the theoretical value of surface pressure 
for close distances of the particular jet as obtained from equation (15). 

The previous experimental data of figure 20 are presented again in figure 21, but 
the jets are considered to have an exit Mach number of 4.79 with an exit half -angle of 
26.5°. From inviscid theory for this Mach number, the two sizes of supersonic nozzle 
would have effective exit diameters of 0.568 in. (1.44 cm) and 1.14 in. (2.88 cm) which 
shifts the experimental data slightly. The equations for a typical exhaust flow case are 
still applicable for the analysis of Roberts and are seen to agree with the experimental 
data as well as before. This slight change in jet-exit Mach number does not significantly 
affect the pressure estimates of the analysis of Roberts. The characteristics plus 
Newtonian pressure estimates, however, change considerably and no longer agree with 
the data. 

The method of characteristics solution for properties in the jet flow field depends 
not only on the jet-exit Mach number but also on the nozzle -exit half -angle. Increasing 
this half -angle from 15° to 26.5° changes the calculated properties in the jet flow field 
considerably and consequently the impingement pressures estimated by using these 
properties. The effective nozzle-exit half-angle for the M g = 5 jet that was calculated 
during the investigation of reference 6 probably resulted from a local effect. Most of 
the gases issuing from the nozzle are affected only by its geometric characteristics. 

The analysis of Roberts neglects this expansion angle as it considers most of the mass 
and momentum of the jet concentrated near the jet axis. 

M g = 7.09 helium jet.- Figure 22 is a summary of the experimental data for both 

sizes of the M g = 7.09 helium jet and includes analytical estimates. Characteristics 
plus Newtonian pressure estimates using the appropriate axial Mach number distribution 
for this jet appear to underestimate the surface stagnation pressures. For the analysis 
of Roberts, the spherical flow equation is applicable. (See fig. 11.) This analysis also 
underestimates the pressures and is seen to coincide with one of the characteristics 
plus Newtonian analyses. For small values of jet height ratio, this analysis should 
underpredict surface stagnation pressures because of the assumption that V = V max 
which is not exactly true near the exit. 

Figure 23 shows the surface stagnation-pressure variation with distance for both 
a stationary and a moving M g = 7.09 helium jet. The small inflection in the data trend 
is possibly a result of some internal -shock phenomenon in the jet plume. The data for 
the moving jet ^0.54 d e /sec^ were taken to determine whether this easier test technique 
might be an acceptable method for obtaining impingement pressures over a range of dis- 
tances. The data for the stationary jet are the same as shown in the previous figure for 
the smaller nozzle. The moving jet appears to provide an expeditious means to obtain 
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impingement pressures over a range of distances provided the continuously decreasing 
ratio of jet -exit pressure to ambient pressure can be overlooked. This pressure ratio 
does not influence the jet exhaust along the jet axis from the jet exit to the intersection 
of the leading characteristic with the axis. For the M = 7.09 helium jet, the leading 

V 

characteristic never intersects the jet axis and thus the axial impingement pressures 
could be determined by using a moving nozzle at any distance as long as the jet remains 
under expanded. The off-axis distribution, however, depends on the pluming of the jet 
which changes with pressure ratio. 

Sonic. jets .- Figure 24 is a summary of the experimental data for both sonic jets 
and includes analytical estimates as well as experimental data of references 7 and 8 for 
a sonic air jet impinging on a flat plate. Pressure estimates from characteristics plus 
Newtonian analyses slightly underestimate the experimental data for the sonic nitrogen 
jet. The analysis of Roberts also underestimates the impingement pressures for this 
jet. This method is not considered to be applicable for low Mach number jets because of 

O 

the assumption that M g » 1 ; however, it was applied to determine its possible useful- 
ness in this range. Two empirical equations which seemed to fair the data are also pre- 
sented in the figure. For the sonic helium jet, the data agree very well with an inverse 
square variation with jet distance. Again, the analysis of Roberts underestimates the 
surface stagnation pressures for this jet. 

Effect of Jet Mach Number 

Figure 25 summarizes the measured surface stagnation-pressure data for the jets 
of the present study. The figure indicates that at a given height ratio, the sonic jets 
cause significantly higher stagnation pressures and hence higher overall impingement 
pressures than supersonic jets with the same gas and jet-chamber pressure. 

Effect of Exhaust Gas Properties 

Figure 26 summarizes the measured surface stagnation-pressure data for the sonic 
jets. The helium jet theoretically produces only about 2.5-percent more thrust than the 
nitrogen jet but experimentally causes about 25-percent higher stagnation pressures for 
the same height ratio and jet-chamber pressure. 

Effect of Surface Curvature 

Pressure distribu tions along cyli nders.- Figures 27 to 30 summarize the distribu- 
tion of impingement pressures along the axis of symmetry of the cylinders compared with 
an analysis for a jet exhaust impinging a flat surface. The comparison of the experi- 
mental data with the analysis of Roberts is made by using the equation 
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( 18 ) 


-^-= (cos 0) K+4 
p stag 

which is the ratio of equation (13) to equation (14). The comparison with a character- 
istics plus Newtonian analysis is made by using the equation 


P S 

P stag 



(19) 


which is the ratio of equation (7) to equation (8). 

Figure 27 is a summary of the impingement pressures for the small M e = 5 
nitrogen jet. The data are plotted against Z/h, which is equivalent to tan 0, without the 
multitude of data points occurring at the jet axis (Z/h = 0). There is some data spread 
for the range of jet distances and model sizes represented; however, the analysis of 
Roberts appears to agree with the data very well. Apparently, the surface curvature of 
the cylindrical models does not appreciably affect the distribution of impingement pres- 
sures along the cylinder axis. The characteristics plus Newtonian analysis is shown for 
two jet heights, a bow shock coincident with the surface being assumed, and does not 
agree particularly well with the experimental data obtained by using the flow character- 
istics of an M g = 4.79, 6 & = 26.5° jet. 

Figure 28 is a summary of the data for the large M e = 5 nitrogen jet. The 
analysis of Roberts agrees with the data nearly as well as that for the smaller jet. 

Figure 29 is a summary of the data for both sizes of the M g = 7.09 helium jet. 
The characteristics plus Newtonian analysis is shown for two jet heights and does agree 
somewhat with the data, but it appears that an empirical equation for which the impinge- 
ment pressure varies as the cos 9 raised to some power would be more applicable. 

Figure 30 is a summary of the data for both sonic jets. Included for comparison 
with the sonic nitrogen jet data are data taken from reference 7 for a sonic air jet 
^d g = 0.50 in. (1.27 cm)j impinging onto a flat surface. These data agree very well and 
are overestimated by the analysis of Roberts. The characteristics plus Newtonian pres- 
sure estimates for the two selected jet heights were identical and slightly underestimate 
the experimental data. The experimental data for the sonic helium jet are also overesti- 
mated by the analysis of Roberts. 

Pressure distributions a round cylinders .- Figures 31 to 33 show pressure distri- 
butions around the cylinders at two jet heights for the different jets together with calcu- 
lations based on the analysis of Roberts by using the equation 
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which is the ratio of equation (10) to equation (14). This equation is easily applied since 
choosing a value of angular position on a given cylinder for a particular jet height uniquely 
determines the angle 9 from geometry as 


Jcos(0 + wjj 


( 20 ) 



0 = tan' 1 

h + R(1 - cos u>) 


( 21 ) 


Figure 31 which shows pressure distributions for the small nitrogen jet and differ- 
ent cylinders compared with the analysis indicates that experimental agreement with the 
analysis depends upon jet height and surface curvature. The analysis of Roberts under- 
estimates the fall-off of impingement pressure with increasing angular position at both 
jet heights for the smallest cylinder, which has the highest surface curvature. For the 
next largest cylinder, the analysis underestimates the pressure fall-off at only the higher 
jet height. The analysis and experimental data agree very well at both jet heights for the 
other cylinders. 

Figure 32 presents similar data for the large M g = 5 nitrogen jet and again indi- 
cates much better agreement between analysis and experiment for the lower jet heights. 
This analysis underestimates the pressure fall-off for the higher jet height for all the 
cylinders. 

Figure 33 shows the distribution of impingement pressures around one of the 
cylinders for the sonic jets at two representative jet heights compared with analysis. 

The characteristics plus Newtonian analysis agrees very well with the data for the nitro- 
gen jet whereas the analysis of Roberts underestimates the fall-off of these pressures. 
The analysis of Roberts also underestimates the pressure fall-off with increasing angular 
position for the helium jet. 

The data from the supersonic helium jet are compared with a characteristics plus 
Newtonian analysis in figure 34 by using the equation 


p s _ p r,l 
p stag p r ,a 



+ “i) 
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( 22 ) 


which is the ratio of equation (6) to equation (8). This analysis indicates some agreement 
although underestimating the fall-off of impingement pressure with increasing angular 
position on the surface. 
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Effect of the Ratio of Jet-Exit Pressure to Ambient Pressure 

The experimental data were obtained at finite values of the ratio of jet-exit pres- 
sure to ambient pressure and it is of interest to determine the effect of this pressure 
ratio. Although the pressure ratios at which the present tests were conducted (see 
fig. 4) were far from being infinite, the jets were not foo far from being fully expanded. 
Figure 35 presents the impingement pressure distribution around one of the cylinders for 
the large = 5 nitrogen jet at a particular value of jet height for three values of 
pressure ratio. The figure indicates that varying the pressure ratio by a factor of 3 did 
not have an appreciable effect on the impingement pressures. The data of reference 8 
are in agreement that significant variations in the ratio of jet-exit pressure to ambient 
pressure at fairly large values of pressure ratio do not influence the magnitude of 
impingement pressures. 


Bow -Shock Characteristics 

Shadowgraphs were taken of the exhaust impingement region of one of the cylinders 
for five different jets from which certain bow-shock characteristics were measured. 

The shadowgraphs were taken in a plane perpendicular to the cylinder axis to determine 
the influence of surface curvature. Representative shadowgraphs are presented in fig- 
ures 36 to 38 for various distances of the jets above the cylinder surface. The figures 
indicate that surface curvature has a decided influence upon the shape of the bow shock. 
When a jet exhaust impinges upon a flat surface, a bowl-shaped bow shock is formed 
which curves away from the surface (see ref. 18); however, the shape of the bow shock 
for a curved surface depends upon the jet, its distance from the surface, and the degree 
of surface curvature. The bow shock is fairly flat when the jets are close to the surface 
and curves toward the surface as the distance between the jets and surface increases 
with the exception of the M g = 7.09 helium jet. The bow shock for this jet has no 
inclination to follow the curvature of the surface as the jet is moved away from the sur- 
face. An explanation for this effect is not apparent; however, this jet is one wherein the 
exhaust leading characteristic never intersects the jet axis. 

The relationship of nozzle size to shock size is shown in figure 39 which is a plot 
of the shock bowl diameters taken from the shadowgraphs of figure 36 for the two sizes 
of M g = 5 nitrogen jet. The data are limited because the shock becomes larger than 
the region being photographed but indicate that the size of the shock is probably directly 
proportional to the nozzle size. An empirical equation which fairs the limited data gives 
the shock size as being equal to twice the distance of the jet from the surface. The size 
of the bow shock should perhaps be expected to increase with increasing ratio of jet-exit 
pressure to ambient pressure since the size of the exhaust plume increases with 
increasing pressure ratio. This condition was shown in reference 18 to exist for an 
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M g = 5 nitrogen jet exhausting onto a flat plate, although the shock characteristics, both 
size and standoff distance, appeared to approach constant values for pressure ratios 
greater than approximately 100. 

The variation of bow -shock standoff distance at the jet axis, as determined from 
the shadowgraphs, is shown in figures 40 to 42. The shock standoff distance for the 
supersonic jets appears to approach a limit with increasing jet distance. The location of 
the bow-shock wave has been suggested as being useful for more closely predicting 
impingement pressures. (See ref. 8.) The bow shock was considered to be coincident 
with the surface for the present study because the usefulness of the analysis should 
depend upon a minimum of information having to be determined. 

CONCLUSIONS 

An experimental study has been made of the impingement of underexpanded jets 
onto various curved surfaces in a vacuum environment and the results have been com- 
pared with available analyses. The study suggests the following conclusions: 

1. The magnitude of axial impingement pressure depends on the exit Mach number 

of the jet and the jet height ratio but not on the shape of the impingement surface. The 
characteristics plus Newtonian analysis adequately predicts this pressure and its varia- 
tion with jet distance from the surface for a range of exhaust jets. The analysis of 
Roberts is easier to apply and predicts this pressure and its variation with jet distance 
from the surface for a range of exhaust jets equally well, provided » 1 (where 

M e is the nozzle -exit Mach number). 

2. A realistic estimate of the axial impingement pressure and its variation with 
jet distance for sonic jets is given by an empirical equation which states that the ratio of 
surface stagnation pressure to nozzle -chamber pressure is equal to the inverse square 
variation of the ratio of distance to jet -exit diameter. 

3. At the same height ratio, sonic jets produce significantly higher axial impinge- 
ment pressures than do supersonic jets with the same gas and jet-chamber pressure. 

4. Sonic helium jets produce slightly higher axial impingement pressures than do 
sonic nitrogen jets at the same height ratio and jet-chamber pressure. 

5. The distribution of impingement pressures depends upon the exit Mach number 
of the jet, the jet height ratio, and the shape of the surface. The characteristics plus 
Newtonian analysis reasonably predicts for a range of exhaust jets the distribution of 
impingement pressures on both planar and curved surfaces. The analysis of Roberts is 
easier to apply and provides about the same degree of accuracy for predicting impinge- 
ment pressures on these surfaces provided M g 2 » 1. 
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I 6. There is no significant difference between the distribution of impingement pres- 

| sures along the longitudinal axis of symmetry of curved surfaces and those along a flat 
plate . 



Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 21, 1966, 
124-08-05-09-23. 
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Figure 2.- Drawing of cylindrical models showing static-pressure orifice locations. All linear dimensions are given in inches and parenthetically in centimeters. 





(c) Convergent-divergent nozzle, d 0 = 1.25 (3.18). 

Section drawing of test nozzles. Linear dimensions are given in inches and parenthetically in centimeters. Because of space 
limitations conversions to the S.l. system of units are presented for only a few representative dimensions. 
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(b) Variation of axial impingement pressure. 
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(c) Variation of the ratio of surface stagnation pressure to 
nozzle -chamber pressure. 
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Figure 5.- Gage response during a typical test. 








Figure 7.- Theoretically calculated contour lines of constant Mach number within nozzle exhaust plume for an M g = 7.09, d Q = 15 jet using helium. 

P c /Poo = 258 X 10 3 . 
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Spherical source flow 



Figure 11.- Flow regimes for conical nozzles delineating the maximum jet-exit Mach number for which the leading characteristic can intersect 

the jet axis. 
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Figure 12.- Distribution of impingement pressures around cylinders for M = 5 nitrogen jet having d = 0.625 in. (1.59 cm). 
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Figure 72- Effect of jet height on surface stagnation pressure for JVL = 7.09, 9 a = 15° jets using helium. 
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Figure 23.- Comparison of axial impingement pressures for stationary and moving = 7.09 jet using helium. 
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Experiment 
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(a) Nitrogen, y = 1.4. 

Figure 24.- Effect of jet height on surface stagnation pressure for sonic jets. 
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(a) Nitrogen, y = 1.4. 

Figure 25.- Effect of jet-exit Mach number on surface stagnation pressure. 
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Figure 26.- Effect of exhaust gas properties on surface stagnation pressure for sonic jets. 
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Nozzle drawn to scale of figure 
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Figure 35.- Effect of ratio of jet-exit pressure to ambient pressure on the distribution of impingement pressures for M g = 5 nitrogen jet. d 

h = 5d ■ d„ = 1.25 in. (3.18 cm). 
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Figure 36.- Concluded, 
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(a) Nitrogen, y = 1.4. L-67-904 

Figure 38.- Shadowgraphs of sonic jets exhausting onto d m = 10-in. (25.4-cm) cylinder. 
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